Wave Propagation Over Parallel Tubular 

Conductors : 

The Alternating Current Resistance 

By SALLIE PERO MEAD 

Synopsis: On the basis of Maxwell's laws and the conditions of con- 
tinuity of electric and magnetic forces at the surfaces of the conductor, the 
fundamental equations are established for the axial electric force and the 
tangential magnetic force in a non-magnetic tubular conductor with parallel 
return. The alternating current resistance per unit length is then derived 
as the mean dissipation per unit length divided by the mean square current. 
The general formula is expressed as the product of the alternating current 
resistance of the conductor with concentric return and a factor, termed 
the "proximity effect correction factor," which formulates the effect of 
the proximity of the parallel return conductor. The auxiliary functions which 
appear in the general formula are each given by the product of the cor- 
responding function for the case of a solid wire and a factor involving the 
variable inner boundary of the conductor. 

In general, the resistance may be calculated from this formula, using 
tables of P :ssel functions. The most important practical cases, however, 
usually involve only the limiting forms of the Bessel functions. Special 
formulae of this kind are given for the case of relatively large conductors, 
with high impressed frequencies, and for thin tubes. A set of curves illus- 
trates the application of the formulae. 

I. Introduction 

WHERE circular conductors of relatively large diameter are 
under consideration, the effect on the alternating current 
resistance of the tubular as distinguished from the solid cylindrical 
form becomes of practical importance. Mr. Herbert B. Dwight has 
worked on a special case of this problem and developed a formula 
for the ratio of alternating to direct current resistance in a circuit 
composed of two parallel tubes when the tubes are thin. 1 As infinite 
sums of infinite series are involved, however, his result is not well 
adapted to computation. 

Mr. John R. Carson has given a complete solution for the alter- 
nating current resistance of two parallel solid wires in his paper 
"Wave Propagation Over Parallel Wires: The Proximity Effect," 
Phil. Mag., April, 1921. The analysis of that paper may readily be 
extended to the more general case of propagation over two tubular 
conductors by a parallel method of development. This is done in 
the present paper. As the underlying theory is identical in the two 
problems, familiarity with the former paper will be assumed and the 
analysis will merely be sketched after the fundamental equations are 
established. 

'"Proximity Effect in Wires and Thin Tubes," Trans. .4. J. E. £., Vol. XLII 
(1923), p. 850. 
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In this paper formulae for the alternating current resistance have 
been worked out in detail with particular reference to the case of 
relatively large conductors at high frequencies and to relatively 
thin tubes. In general the auxiliary functions involved are expressed 
as the product of the corresponding functions for solid wires by a 
correction factor which formulates the greater generality due to 
the variable inner boundary of the conductors. As far as possible 
the symbols are the same as in the solid wire case but refer now to 
the system of tubular conductors. Primes are added where the 
letters denote the corresponding functions for the solid wire case. 
This will hardly lead to confusion with the primes used in connection 
with the Bessel functions to denote differentiation. 

The general solution is developed in section II. The alternating 
current resistance of one of the tubular conductors is expressed as the 
product of the alternating current resistance of the conductor with 
concentric return and a factor which formulates the effect of the 
proximity of the parallel return conductor. Section III is a sum- 
mary of the general formula, special asymptotic forms and forms 
for thin conductors. 

II. Mathematical Analysis and Derivation of Formulae 

We require the expression for the axial electric force, E z , in the 
conductors. Since the tubular conductor does not extend to r = 0, 
the electric force must be expressed by the more general Fourier- 
Bessel expansion, 

eo 

E t = ^rA„[J„(p)+\„K„( P )] cos nd, 



n=0 



where 



= £ = xiy/i when r=a 

= f = yi v * when r = a, 

a and a being the outer and inner radii, respectively, of the con- 
ductors. The additional set of constants X , Xi . . . X» is to be deter- 
mined by the conditions of continuity at the inner boundary of the 
conductor. It is necessary to satisfy the boundary conditions at the 
surface of one conductor only, since the symmetry of the system 
insures that they will then be satisfied at the surface of the other also. 
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In the dielectric space inside the tube where r<a, the axial electric 
force may be written 

E z = V C„J H (p) cos nd, (1) 

n = 

or replacing the Bessel functions by their values for vanishingly 
small arguments, 

oo 

E s =^D n r n cos nd (2) 

«=o 

where D , D\ . . . D„ are constants determined by the boundary condi- 
tions. Applying Maxwell's law relating the normal and tangential 
magnetic forces H r and H e to the axial electric force, gives 

»ia>H g = ^ ]T A„ [Jn'( P ) + *nK'„( P )] cos nO, (3) 

n=0 

fiiuHr = - ^ A n [J„(p) +\JC„(p)] sin w0, (4) 

n = 

for the space inside the conductor, and 

iwH g = y ^ nD„r"- 1 cos nO, (5) 



n=0 



IwHr = 2 nr) » r "~ l sm ne < (6) 



*=o 



for the inner dielectric (m = 1)- Equating the two expressions for 
the tangential magnetic force H e and for the normal magnetic in- 
duction fxHr term by term at the surface r=a, 

[f/» / (f)-Mw/«(f)]+x„[r^» / (f)-M«^n(rt]=o. (7) 

Whence, for the practically important case of non-magnetic con- 
ductors in which fi = l, we have 



and 



\ -Ai-n(f) /on 



Es = 2 A * [ JM ~ I^M KM ~1 cos "* (9) 



« = 
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In the subsequent analysis /„ (£) of the solution for the solid wire 
case is replaced by 

M& - j^fy^) = M »W> do) 

and J n ' (£) is replaced by 

Jn'(k)-^^jKn'(0=M„'(k). (11) 

Otherwise the formulation of the alternating current resistance of 
the conductor proceeds exactly as in the solid wire case. For the 
electric force at the surface r=a in the conductor, we write 

Ei^AolMoiD+hiM^) cos-0+Mf«(f) cos 20+ . . .] (12) 

and determine the fundamental coefficient A in terms of the current 
in the conductor. The resistance R of the tubular conductor per 
unit length is defined as the mean dissipation per unit length divided 
by the mean square current where the mean dissipation is calculated 
by Poynting's theorem. Accordingly, we get 

2w*(M e (i) , 1 Viz n^-^mm Mn '^l run 
^ Rea V TO + 2l |/4 TO con, 'Il ' (13) 

n=l ' 

To determine the harmonic coefficients hi . . . h n or Ai . . . A„, 
the total tangential magnetic force and the total normal magnetic 
induction at the outer surface of a conductor are expressed in terms 
of the coordinates of that conductor alone, and the conditions of 
continuity at the surface are applied. This leads to the set of equations 

qn = (-l)*2p n k»--f0£ [P nk"^ n (q) (14) 

w=l,2,3 ... =o 

where 

V /n w! z, («+l)!,, . 

<r n = ^Mn {k)-nnM n (?))/? M 7 ^ 

Pn = (*!£*' (Q - nfxM n (0)/(M n '^) +»/«Jf.(0) l 

([n = (Tnhn, 

C 
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When the permeability is unity, the solution, to the same order of 
approximation as in the solid wire case, is 

l+XJCi(D//i(0 



w-* ,+ * 



i+x„A'„_iU)//„-.a) 



! Mi+2»gfe 2 A"-') ( 16 ) 



where 

_ y / 2~ p[«i(^+^)~ z, iO < '>~ ; '°)l~gt"i( ?/o ~ t ' )+ t ' 1 ^ + t ' ^ (17) 

S_ x th 2 +v 2 

p n = (-l) n 2£ B 5 n , w = l,2 . . . po, 

1-VT-(W 

5_Z (2fc) 2 
Since the resistance R of an isolated tubular conductor is given by 

i^Real—^-^ (19) 

equation (13) becomes equation (I) of the formulae in the next section. 
This is the general solution for the case of non-magnetic conductors. 
In general R may be calculated from this formula and tables of 
Bessel functions. The ber, bei, ker and kei functions 2 and the recur- 
rence formulae are sufficient to evaluate the Bessel functions but 
the process is long. In the most important practical cases, the 
conductors are rather large and the applied frequencies fairly high. 
When this is true as well as when the tubes are very thin the formulae 
usually involve only the limiting forms of the Bessel functions. These 
special results are given in the next section. 

III. Alternating Current Resistance Formulae for 
Non-Magnetic Conductors 

The symbols used are : 

a=outer radius of conductor in centimeters, 

a = inner radius of conductor in centimeters, 

r = interaxial separation between conductors in centimeters, 

k=a/c 

X = conductivity of conductor in electromagnetic c.g.s. units, 

2 A convenient table of these functions for arguments from to 10 at intervals of 
1 is incorporated in Mr. Dwight's paper "A Precise Method of Calculation of 
Skin Effect in Isolated Tubes," J.A.I. E. E., Aug., 1923. 
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H = permeability of conductor in electromagnetic c.g.s. units, 
u = 2tt times frequency in cycles per second, 

x =a\/4ir\oj 
y = ay/Airhw 
£ = xi\/l 
£ = yi\/i 

x«=--/„+i(r)/A'„ + ,a) 

= Bessel function of first kind of order n and argument xiy/I, 

dJn(Z) 



Jn(£) = 
Un'+iv„' = 



dJ„tt) 



dx 
K (£) = Bessel function of second kind of order n and argument xiy/T, 

R = resistance per unit length of tubular conductor with parallel 
return, 

R = resistance per unit length of tubular conductor with con- 
centric return in electromagnetic c.g.s. units, 

C = proximity effect correction factor, 

R = CR Q . (I) 

The auxiliary functions involved are: 

ID d / ( '1 ntioHo -\-v v ' \ 

8 Ro = Rom I 1 t — ) (20^ 

\ mUoVo —u o v l 

where 

ft ,.iEaa^ ( 2i) 

a \ 7rX Ui 2 +v x - ' 

= resistance of solid wire with concentric return, 

i±sMMM£ (22) 

a - o'-h S 1 - g [«l(«o — Po)+Pl(«o+g'o)l I , 9 „v 

g gP \ l phh(u +v )-v 1 (u o -v ))\' W 

1 The ratio R„/R'„ oscillates about unit}' which it approaches more and more 
closely as the frequency increases. It is due to the fact that the phase of the current 
in the inner portion of the solid conductor may be such as to oppose the current 
in the outer portion, that the resistance of the solid conductor may be greater than 
that of the tube even though the heating effect in the latter is the greater. 
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where 



, _ \/2 H i ( »o + Vo) — Ui( II o — Vo) 
* ~ X U 2 -\-Vo 2 



l+Xi£itt)M($) 

b„ UnU,/+V„V n '' 



W 

where 



, UnVn — 1(n V„ 

W„ = 



where 



n-l 

For large values of the argument 



(24) 



^ +, « = l + X 1 Xo(€)//o(«)' (25) 



"" W " | l + X„K„-iU)M<-i(£) | 2 \ O- UnVn'-Un'vJ' K * 

(27) 



o, + f6,-(l + X,^>)con j .(l + X,^) > (28) 

1-Vl~(2fe)» ( 29) 

(2fc) 2 
The formula for the correction factor C is then 

aK y tt\ 

00 

5i=Vw«ife 2 " 5 2 », (30) 

M=l 
OO 

.5 2 =^my„P"5''+ 1 . (31) 



&-*'[«•— (l- -^ )] (32) 



V2.V 
and the correction factor is 

C=1+2 _^( 5l -^(l4)>'5,) (III) 

When .v and y are both large quantities, the auxiliary functions are 
as follows, provided terms of the second order in l/x and \/y are 
negligible, n in d and h below being equal to the number of terms 
in which Si and S 2 converge to a required order of approximation. 

With the notation 

cos = cos v2 (x — y) , 
sin = sinV2(#— y), 
exp = exp [ — s/2(x —y)] t 

, l + [(l+q) sin-(l-g) cos] exp -a exp 2 . . 

Ko-Ko 1 _[( 1 _ 6 ) sin + (i + 6 ) cos ] exp _ 1 _6exp 2 K ° ' 
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where 

a = l — 



6 = 14- 



2-s/2x- 2V2.v' 
3 3 



aR 



where 



2\/2x 2\/2y 



_ , ! + [(! — c) cos— (1 + c) sin] exp — c exp 2 
~ g l-[(l + c)cos+(l-c)sin]exp+cexp 2 ' 



(34) 
(35) 



1 15 

c = l — 



where 



2\/2.t: 2\/2y 

g'=--v/2>, (36) 

-., ^ ?f . . l-[(l-^)cos-(l+^)sin]exp-rfexp 2 

Wn Wn l-[(l-f-//.)cos4-(l-/0sin]exp4-Aexp 2 ' Cd/j 



, , . 4w 2 -l 4(w + l) 2 -l 
a = 1 



fc = 14 



2\/2* 2y/2y ' 

4(n-l) 2 -l 4(w4-l) 2 -l 



2y/2x 2\/2y ' 

Wn= V5~-*r- (38) 

At frequencies sufficiently high to afford practically skin conduction, 
the following formulae indicate the way in which the resistance of 
the tubular conductor approaches its limit, the resistance of the 
solid wire. 

R, = R.' \ + l sine * P , (39) 

1 — 2 cos exp 






_1 I U> 

aR 



C=C m (l-A/x), (IV) 

Cm_ l^P?' (40) 

A-n^/n- W Jt ■ o» (l-^ ! 1-2 sin exp j 

^- 2 V2 r^v| 1 + 2 ^ 2 (I3^)T 1 _2cosex P r (41) 
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When the conductors are very thin tubes, i.e., thin as compared to 
the radius, (a-a)'a is necessarily small and, in general, x-y is 
small. Of course, when the frequency is high enough, x-y becomes 
large in any case. When this is true with respect to thin tubes, how- 
ever, x and y will usually be large enough to make the asymptotic for- 
mulae applicable; but, if x-y is small, the approximations 

/.tf)-/.(0-tt-r)/-'(D+^p/."(D. 

if»(f)^»(a-(?-s)ivV(a,+ tfW'(£), 



reduce the correction factor to 



£-*■ 



C-l+2flf| JP W'^ - 2fe 2 ** F V^"5»+»n^- } (V 

n=l " "=l 



where (3 = 



a — a 



f = 



(l+/3/2) 2 _c 2 



1+/3+/9 2 do' 



A, = /3V,, 2 -r-^W, 



(»+l)(n+2) 



4,-l+(*+l)04 

and the resistance with concentric return to 

1 1+/3+/3 2 



/3 2 . 



*o = 



27r>a(a-a) 1+/3/2 ' 



(42) 



l/2irXa(a — a) is, of course, the direct current resistance of a very thin 
conductor. 

If (a — a)/a is very small and negligible compared with 2n/x~, where 
n is the number of terms in which the series of (V) converge to a 
required order of approximation, 



c-i+fC-TT 



(i-^) \X 



^-+2£ 2 slog(l-fc 2 s) 



-«jlog (1-^ + 2^ j 



(VI) 
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As a check on formulae (V) and (VI), the limiting cases may be 
arrived at directly as follows. If the conductors are thin tubes, 
the harmonic coefficients are given by 

//„=:(- l)»+>2&" 



When £ is very large 

-(-1)»2*«j«, (44) 



and 



M Mo 



so that 

i «« r 

l+£ 2 s 2 



c-w[i + *2i*.p£:«*g] 



1-w 



(46) 



the same result as for the corresponding limiting case of a solid con- 
ductor. 

On the other hand, if £ is not large and £ — f is very small, 

hn=(-iy+^^-n, (47) 

K- 1 ' (48) 



m' 



7>7 



(49) 



so that 
and 



M fl ' *(*-y)' 

C = l, (50) 

R = Ro=R d ,., (51) 
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where R dx . is the direct current resistance of the thin tubular con- 
ductor. Eqs. (46) and (50) agree with the corresponding limits of 
formulae V and VI respectively. 

The curves of the accompanying figure do not pretend to represent 
the proximity effect correction factor with precision. They are, how- 
ever, accurate for thin tubes, and indicate the order of magnitude 
of the factor for various values of the thickness of the tubular con- 
ductor and show the nature of its variation with respect to the applied 
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frequency. They are computed from formula (V) which is valid 
for quite high frequencies when the tubes are thin. When the thick- 
ness of the tubes is greater, however, the range of validity with respect 
to frequency is smaller, the dotted portions indicating a doubtful 
degree of precision. It was previously pointed out in connection 
with formula (IV) and is immediately deducible from physical con- 
siderations, that all of the curves eventually coincide with the curve 
for the solid wire which approaches the value 1.155 asymptotically. 
As a simple application, suppose the resistance is required of a 
tubular conductor with an outer radius of 0.4125 cm. (that of No. O 
gauge A.W.G. copper wire) whose resistivity is 1696.5 electromagnetic 



338 BELL SYSTEM TECHNICAL JOURNAL 

units per cm., where there is an equal parallel return so situated that 
k = 0.25 and a frequency of 5,000 cycles per second is applied to the 
circuit. Then m = \/4ir'Ko) = 15.26 and x=wa = 15.26X0.4125 = 6.30. 
When the ratio of the thickness of the conductor to the radius is 
greater than about 0.01 the proximity effect correction factor C is 
appreciable. If the ratio is 0.05, reading C from the curves, gives 
C= 1.064. From formula (42), i? = 5.24 ohms per mi. which makes 
the resistance i? = 5.53 ohms per mi. 



